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Abstract 

Recent higher-order results for precision observables in the Standard Model 
(SM) and the Minimal Supersymmetric Standard Model (MSSM) and for the 
neutral CP-even Higgs-boson masses of the MSSM are summarized. Further- 
more a brief discussion of the technical aspects of evaluating higher-order cor- 
rections in the electroweak theory is given. In the SM, results for the Higgs- 
mass dependence of the precision observables are analyzed. The exact two- 
loop results for the Higgs-mass dependence of the fcrmionic contributions are 
compared with the results of an expansion in the top-quark mass up to next- 
to-leading order. In the MSSM, results for the leading two-loop contributions 
to the precision observables and to the masses of the neutral CP-even Higgs 
bosons are discussed. The latter are compared with results obtained by renor- 
malization group calculations. 
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Recent higher-order results for precision observables in the Standard Model (SM) 
and the Minimal Supersymmetric Standard Model (MSSM) and for the neutral 
CP-even Higgs-boson masses of the MSSM are summarized. Furthermore a brief 
discussion of the technical aspects of evaluating higher-order corrections in the 
electroweak theory is given. In the SM, results for the Higgs-mass dependence of 
the precision observables are analyzed. The exact two-loop results for the Higgs- 
mass dependence of the fermionic contributions are compared with the results of 
an expansion in the top-quark mass up to next-to-leading order. In the MSSM, 
results for the leading two-loop contributions to the precision observables and to 
the masses of the neutral CP-even Higgs bosons are discussed. The latter are 
compared with results obtained by renormalization group calculations. 

1 Introduction 

By comparing the electroweak Standard Model (SM) and its extensions, most 
notably theJVlinimal Supersymmetric Standard Model (MSSM), with the pre- 
cision data El it is possible to test the theory at its quantum level, where all 
parameters of the model enter. In this way one is able within the SM to infer 
indirect constraints on the mass of the Higgs boson, which is the last missing 
ingredient of the SM. Within the MSSM the comparison with the precision 
data allows to constrain mainly the parameters of the scalar quark sector. 
Eventually the precision tests of the electroweak theory could hint towards a 
distinction between the SM and Supersymmetry via their respective virtual 
effects. At present the bounds that can be obtained from this analysis on the 
parameters of the SM and the MSSM are still relatively weakld. In order to 
improve this situation, a very high accuracy both of the measurements and the 
theoretical predictions is needed. 

While QCD corrections to Ar and the Z-boson observables sin 2 9 e g and Ti 
are known within the SM at C(aa s )EI and 0(aa 2 )tj, only partial information 
is available about the electroweak two-loop contributions. Beyond one-loop 
order the resummations of the leading one- loop contributions are known ltd, 
the leading and next-to- leading term in an expansion for. asymptotically large 
values of the top-quark mass, m t , have been evaluatedEffl'ta, and also the lead- 
ing term of an asymptotic expansion in the Higgs-boson mass, Mh, has been 
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derived □. The terms obtained in the mt-expansion were found to be numer- 
ically sizable, and the next-to-leading term turned out to be about equally 
large as the (formally) leading termQ. Exact results have been derived for the 
Higgs-mass jienmdence of the fermionic two- loop corrections to the precision 
observablesBlU. 

In order to treat the MSSM at the same level of accuracy as the SM, 
the one- loop results for Ar and the Z-boson observables LL3 have to be sup- 
plemented by higher-order contributions. Recently—the QCD corrections to 
the p parameter in the MSSM have been evaluated Ei3, which incorporates the 
leading two-loop contribution to the precision observables. 

Besides the indirect tests of the MSSM via the precision observables, there 
exists a very stringent direct test of the model, since it predicts the existence 
of a light Higgs boson, which at the tree level is restricted to be lighter than 
the Z boson. A precise prediction for the mass of the lightest Higgs boson, 
mh, in terms of the relevant SUSY parameters is crucial in order to determine 
the discovery and exclusion potential of LEP2 and the upgraded Tevatron. 
If the Higgs boson exists, it will be accessible at the LHC and future linear 
colliders, where then a high-precision measurement of the mass of this particle 
will become feasible. A precise knowledge of the mass of the heavier CP-even 
Higgs boson, toh, is important for resolving the mass splitting between the 
CP-even and -odd Higgs-boson masses. 

The tree-level bound on mh is strongly affected by the inclusion of radiative 
corrections, which yield an .ni^pe, 
yond the one-loop resultsl 14 ' 1 ' 1 



ound of about 130 GeV. In order to go be- 
, renormalization group (RG) methods 
been applied for obtaining leading logarithmic higher-order contributionsEi 
In the effective potential approach the leading QCD corrections have been cal- 
culatedo. Up to now phenomenological analyses have been based either on the 
results of the RG improved one-lopppeffective potential approach or on the 
complete one-loop on-shell results liaO, which differ by up to 20 GeV in mh- 
Recently a Feynman-diagrammatic calculation of the leading QCD corrections 
to the masses of the neutral CP-even Higgs bosons has been performed Eilua. 
The result obtained in this way contains thedeading two-loop corrections, the 
full diagrammatic one-loop on-shell result 113, and further improvements tak- 
ing into account leading electroweak two-loop and leading QCD corrections 
beyond 0(aa s ). 

In this paper some recent higher-order results in the SM and the MSSM 
are summarized. Before discussing these results in some detail, a brief overview 
is given over technical issues involved in the calculation of two-loop corrections 
in the electroweak theory. 
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2 Techniques for calculating higher-order corrections in the Stan- 
dard Model and the MSSM 

The evaluation of higher-order corrections in the electroweak theory to the pre- 
cision observables and to the predictions for the MSSM Higgs-boson masses 
involves as a main technical obstacle the evaluation of two-loop 2-point func- 
tions. In order to express the results for the precision observables in terms of 
the physical masses A% and Mz of the gauge bosons, the respective 2-point 
functions have to be evaluated on-shell, i.e. at non-zero momentum transfer, 
while vertex and box contributions in processes with light external particles 
can often be reduced to vacuum integrals. As a consequence of the many dif- 
ferent mass scales present in the electroweak theory, the evaluation of diagrams 
at the two-loop level is in general very complicated. 

Besides the large number of contributing Feynman diagrams, problems 
encountered in such a calculation are due to the complicated tensor structure 
of the diagrams and to the fact that the scalar two-loop integrals are in general 
not expressible in terms of polylogarithmic functionscij. Further complications 
are related to the need for an adequate regularization and for a renormalization 
at the two-loop level, which has not yet been worked out in full detail. 

Concerning the regularization, two schemes are commonly applied in cal- 
culations_within the electroweak theory, namely Dimensional Regularization 
(DREG) B@ and Dimensional Reduction(DRED) @. In DREG the regular- 
ization is performed by analytically continuing the space-time dimension from 
4 to D. This prescription preserves the Lorentz and the gauge invariance of the 
theory, apart from problems related to the treatment of 75 in dimensions other 
than 4. In Supersymmetric theories, however, a D-dimensional treatment of 
vector fields leads to a mismatch between the fermionic and bosonic degrees 
of freedom, which gives rise to a breaking of the Supersymmetric relations. 
This led to the development of DRED. In this scheme only the momentum 
integrals are treated D-dimensional, while the fields and the Dirac algebra are 
kept 4-dimensional. DRED involves potential ambiguities related to the treat- 
ment of 75 , and its application to non=supersymmetric theories turns out to be 
problematic (for a recent review see E3) . In a naive application of both regu- 
larization schemes (without an appropriate shift in the couplings which relates 
the two schemes to each other) at one-loop order differences in the finite parts 
of the Feynman diagrams (proportional to {D — 4)°) arise, while at two-loop 
order both the finite parts and the divergent part proportional to (D — 4)" 1 
are different. 

As already mentioned, a problem closely related to the one of finding an 
adequate regularization for the considered calculation is the treatment of 75 
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in D dimensions. The consistent prescription according to E§E3 leads to a 
breaking of the Ward identities. As in all cases where one works with a non- 
invariant regulator, this breaking has to be compensated by the introduction 
of additional counterterms that restore the symmetries. In order to investi- 
gate the invariance properties of the regulators in Supersymmetry, a detailed 
study of the supersymmetric Ward identities is necessary. A new approach 
that could be helpful for applications in the electroweak theory is differential 
renormalization E3. 

The results presented in this paper have been obtained using several com- 
puter-algebraic tools E3. The generation of the diagrams and counterterm 
contributions has been done with the help of the computer-algebra program 
FeynArtsM. For the calculations in the MSSM, the relevant part of the MSSM 
Lagrangian has been implemented into FeynArts. The program inserts propa- 
gators and vertices into the graphs in all possible ways and creates the ampli- 
tudes including all symmetry factors. The calculation of the two- loop diagrams 
and counterterms was performed with the package TwoCalcS^. It is based on 
a general algorithm for the tensor reduction of two-loop 2-point functions, and 
reduces the amplitudes to a minimal set of standard scalar integrals. Prop- 
erties like the gauge-parameter dependence or the validity of Ward identities 
can directly be read off from the algebraic result obtained with TwoCalc. 

For the calculations in the SM, DREG has been applied, while the re- 
sults within the MSSM presented below have been obtained using DRED. For 
the calculations described in this paper the use of an anticommuting 75 in 
D dimensions was possible without encountering inconsistencies. The renor- 
malization has been carried out using the on-shell scheme. Within the SM 
a two-loop renormalization was necessary in the gauge-boson sector. In the 
MSSM a two-loop renormalization in the Higgs-boson sector and a one-loop 
renormalization of the gauge-boson and scalar quark sector had to be per- 
formed. Those two-loop scalar integrals for which no analytic expression in 
terms of polylogarithmic functions can be derived have been evaluated numer- 
ically with one-dimensional integral representationsEj. These allow a very fast 
calculation of the integrals with high precision without any approximation in 
the masses. 

3 Higgs-mass dependence of precision observables in the Standard 
Model at two-loop order 

In order to study the Higgs-mass dependence of the precision observables Ar, 
sin 2 9 c g and T; at two-loop order, we consider subtracted quantities of the form 

a su btr(Aftj) = a(M H ) - a(M&), where a = Ar, sin 2 9 cS , Ti, (1) 
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which indicate the shift in the precision observables caused by varying the 
Higgs-boson mass between M H and Mr- In the analysis below the Higgs- 
boson mass is varied in the interval 65 GeV < Mh < 1 TeV. The quantity Ar, 
which determines the relation between the vector-boson masses in terms of the 
Fermi constant, is derived from muon decay Ej. It is defined according to 

The leptonic effective weak mixing angle and the leptonic width of the Z boson 
are determined from the effective couplings of the neutral current at the Z- 
boson resonance, J^ c = (v^G^Ml) ^ 2 5y7 M — g^l^ls , according to 

(3) 

Potentially large Mn-dependent contributions are the ones associated with 
the top quark due to the large Yukawa coupling ttH, and the corrections pro- 
portional to Act. Further contributions are the ones of the light fermions 
(except the corrections already contained in Act) , and xuirely bosonic contri- 
butions. The latter are expected to be relatively smallliJ. We therefore have 
focussed on the Higgs-dependent fermionic contributions to the precision ob- 
servables, for which we have obtained exact two-loop results EjO. 

Fig. [l] shows the Higgs-mass dependence of the two-loop corrections to Ar 
associated with the top/bottom doublet, with An, and with the light fermions. 
The dotted line furthermore indicates the Higgs-mass dependence of the lead- 
ing TOj-termcl in the top/bottom contribution. The two-loop top/bottom con- 
tribution gives rise to a shift in the W-boson mass of AM^ P subtr / 2 » (Mh = 
1000 GeV) 16 MeV, which is about 10% of the one-loop contribution. The 
Higgs-mass dependence of this contribution turns out to be very poorly approx- 
imated by the leading m^-term; the contribution of the latter even enters with 
a different sign. It can be seen from Fig. [j] that the two-loop contributions to 
a large extent cancel each other. The contribution of the light fermions yields 
a shift in M\y of up to 4 MeV. In total the two-loop contributions lead to a 
slight increase in the sensitivity of Ar to the Higgs-boson mass compared to 
the one-loop case. Similarly as for Ar, also for the Higgs-mass dependence of 
sin" 9 c g and T/ large cancellations occur between the two- loop contributions. 
For these observables the higher-order contributions decrease the sensitivity to 
the Higgs-boson massEJ. The Higgs-mass dependence of Myv, sin 2 9 e g and T; 
is shown in Fig. || together with the experimental values of the observables. 
The theoretical values for M H = 65 GeV are taken fromn. 
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Figure 1: Higgs-mass dependent fermionic contributions to Ar at two-loop order. The 
different curves show the contribution from the diagrams involving the top/bottom doublet 
(Ar t b), the contribution proportional to An (ArAo), the contribution of the light fermions 
(Arjf ), and the approximation of the top/bottom correction by the leading term proportional 
to m\ (Ar(m*)). 

In Tab. | the Higg s-mass dependence of Myy and sin 2 9 e s based on the 
exact result for the fermionic contributions is compared with ±he results of the 
expansion in the top-quark mass up to 0(G 2 m 2 Af|) given inQ (with the input 
parameters as in H) .^| Over the range of Higgs-mass values from 65 GeV to 
1 TeV the difference between the results amounts to about 4 MeV for M\v and 
to about 7- 10~ 5 for sin 2 c g. The difference in the prediction for sin 2 9 e g turns 
out to be mainly induced by the difference in the prediction for Mw(Mh). 
Evaluating sin 2 e ff, S ubtr(M H = 1 TeV) using the value for M W (M H = 1 TeV) 
fromH instead of our result (see Tab. wields a value for sin 2 # e ff,subtr(MH = 
ITeV) that differs from the one given inla only by about 2 • 10" 5 . 

In the comparison of Tab. |l| besides the difference between the exact treat- 
ment of the top-quark contributions and the expansion up to next-to-leading 
order in m t further effects enter. In lI the Higgs-mass dependence of light- 
fermion contributions only arises from reducible contributions that are gen- 
erated by Dyson summation, while in our results above also the Higgs-mass 
dependence of the light-fermion contributions is taken into account exactly. 

b It should be noted that in contrast to the numbers given in El the contribution of the light 
fermions is included here. 
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Figure 2: The Higgs-mass dependence of the prediction for Mw, sin 2 9 e g and T; is shown 
for different values of nit. The experimental values of the observables are also indicated. 
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Table 1: The Higgs-mass dependence of J% and sin 2 c g based on the exact result for the 
fermionic contribution (left column) and on the result of the expansion in mt (right column). 
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Further differences are due_tp a different treatment of higher-order QCD and 
electroweak corrections. InEJ the effect solely caused by the difference between 
the exact treatment of the top-quark contributions and the expansion in mt 
has been analyzed in detail, and it has been found that over the full range of 
Higgs-boson masses this effect alone amounts to a difference of about 2 MeV 
in Myy and 3 ■ 10~ 5 in sin 2 9 e f[. In the case of My? the three other sources 
of differences mentioned above, namely the light-fermion contribution and the 
treatment of higher-order QCD and electroweak corrections, all individually 
give rise to a difference of 2 — 3 MeV over the full range of Higgs-boson masses, 
which combine with the effect of the m t -expansion to the difference of 4 MeV 
listed in the last entry of the first table in Tab. |l|. 

4 QCD corrections to precision observables in the MSSM 

The leading two-loop corrections to the electroweak precision observables in 
the MSSM enter via the quantity Ap, which is given in terms of the transverse 
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parts of the two-loop Z-boson and W-boson self-energies as 

£ z (0) £ w (0) 
Ap= lrf 2 AP ' (4) 

The contributions to Ap of 0(o!a s ) have been evaluated in El. They can 
be separated into the contribution of diagrams with gluon exchange, which 
dominates in general (it is of the order of 10-15% of the. one-loop result) and 
for which a very compact analytic result can be derivedlia, and into the gluino- 
cxchange contribution, which goes to zero for large gluino mass. 

The results oflifj have recently been extended by including also the effects of 
mixing in the scalar bottom sector, which previously had been neglected, and 
by evaluating the gluon-exchange contributions to the precision observables 
Ar, sin 2 # e ff and Ti, thus going beyond the Ap approximation In Fig. ||the 
gluino-cxchangc contribution is shown for tan j3 — 40 and large mixing in the 
scalar bottom sector as a function of the common scalar mass parameter 

m s = Mi = Mi = Mr = Mr . (5) 

The Mq L/R are the soft SUSY breaking parameters in the diagonal entries 
of the mass matrices of the scalar top and bottom quarks. The off-diagonal 
entries read m q M^ R with M^ R = A t - fi cot (5 , M^ R = A b - [it&nfl , seeEI. 
The full result is compared in Fig. || to the case where mixing in the scalar 
bottom sector is neglected (nib = 0). The figure shows that effects of b mixing 
can be sizable for large tan/3 and small gluino masses. 

The leading contribution to Ar in the MSSM can be approximated by the 
contribution to Ap according to Ar w —c^/s^Ap, where = 1 — = 
M-^/M|. The gluon-exchange contribution to Ar of a squark doublet is 
given by 



Ar susY _ W(Q) 4 (SMl SM^\ S w (0)-^ 
Ar gluon - 11 (0) - — y-j^- -^)+ M~w ' 



where SM^ = ReE w (M^), 5M| = ReS z (Af|), and IF, E w , and S z de- 
note the transverse parts of the two-loop gluon-exchange contributions to the 
photon vacuum polarization and the W-boson and Z-boson self-energies, re- 
spectively, which all are understood to contain the subloop renormalization. 
This result is shown together with the Ap approximation in Fig. ||. The two- 
loop contribution leads to a shift in the W-boson mass of up to 20 MeV for 
low values of nig in the no-mixing case. If the parameter M^ R is made large 
or the relation eq. <M) is relaxed, much bigger effects are possible Eil As can 
be seen in Fig. ^, the Ap contribution approximates the full result rather well. 
The two results agree within 10-15%. 
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Figure 3: Contribution of the gluino-exchange diagrams to Ap for tan (3 = 40 and large 
mixing in the scalar bottom sector. The full result is compared to the case where mixing in 
the scalar bottom sector is neglected (m^ = 0). 
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Figure 4: Contribution of the gluon-exchange diagrams to Ar SUSY as a function of the 
common scalar mass parameter m,q for the scenario of no mixing (M^ R = 0) and large mixing 
(M^ R = 200 GeV) in the t sector. The exact result is compared with the approximation 
derived from the contribution of Ap. 
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Figure 5: The mass of the lightest Higgs boson in the MSSM in terms of the physical 
parameters , and where Am ( - = — mj . The scenarios Am,- = GeV, 0j = 
(no mixing) and Am t - = 340 GeV, 6f = — 7r/4 (maximal mixing) are shown. 



5 Diagrammatic two-loop results for the masses of the neutral el- 
even Higgs bosons in the MSSM 

In0@ a Feynman-diagrammatic calculation of the leading QCD corrections to 
the masses of the neutral CP-even Higgs bosons has been performed. The result 
obtained in this way contains the leading two-loop corrections, the complete 
diagrammatic one-loop on-shell resultEj, and further improvements taking into 
account leading electroweak two-loop and leading QCD corrections beyond 
0(aa s ). The calculation of rah and mu has been carried out for arbitrary 
values of the parameters in the Higgs and scalar top sector of the MSSM. The 
results have been implemented into the Fortran program FeynHiggszH. 

In Fig. IH the mass of the lightest Higgs boson is shown as a function of 
the mass of the heavier scalar top quark, , for different values of mj and 
the mixing angle in the scalar top sector, % The tree-level, the one-loop 
and the two-loop results for mi, are given for the case of no mixing (0j = 0) 
and degenerate t masses (Amj = m t ~ 2 — roj = 0), as well as for maximal 
mixing (#j = — 7r/4) and a mass difference Amj « 340 GeV that gives rise to 
the largest values of mk_/the corresponding values for Ap SUSY lie within the 
experimentally favoredES region Ap SUSY < 1.3 • 10~ 3 ). As can be seen from 
the figure, in both cases the two-loop corrections give rise to a large reduction 
of the one- loop on-shell result of up to 15 GeV. 
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Figure 6: The maximally possible value for m h as a function of tan/3 for rriq = 1000 GeV 
and three values of mt. 

In order to derive an upper bound for mh within the MSSM as a function 
of tan j3, we have performed a scan over the other parameters such that the 
corresponding value of mh is maximized. Fig. ^ shows the maximally possible 
value for mh in the region tan/3 < 5 for rriq = 1000 GeV and three values 
of the top-quark mass, m t = 173.8 GeV, 178.8 GeV, 183.8 GeV (the current 
experimental value of mt and values one and two standard deviations above 
it). It can be seen that the maximal values of mh in the particularly interesting 
region of tan /3 < 2 are at the edge of the LEP2 reach, which is expected to be 
roughly 105 GeV. The question whether a full coverage of this region is possible 
depends sensitively on the range of m t values employed for evaluating m^. 

In Fig. the results of the Feynman-diagrammatic calculation for mh are 
compared with the results obtained by RG methodslij. The mass of the light- 
est Higgs boson is shown as a function of M t / rriq for the two scenarios with 
tan/? =1.6 and tan/3 = 40. Good agreement is found for vanishing mixing 
in the scalar top sector [M^ R = 0). Sizable deviations occur, however, when 
mixing in the t sector is taken into account. They reach about 5 GeV for mod- 
erate mixing and become very large for \M^ R /rriq\ > 2.5. The maximal value 
for mh in the diagrammatic approach is reached for M^ R /rriq w ±2, whereas 
the RG results have a maximum at M t /rriq w ±2.4, which corresponds to the 
maximum in the one-loop result. In the case of positive M^ R , the maximal 
values for mh reached in the diagrammatic calculation are up to 5 (3) GeV 
larger than the ones of the RG method for tan/3 =1.6 (40). A comparison of 
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Figure 7: Comparison between the Feynman-diagrammatic calculation and the results 
obtained by renormalization group methods. The mass of the lightest Higgs boson is shown 
for the two scenarios with tan/3 =1.6 and tan/3 = 40 as a function of /m,q. 



the diagrammatic and the RG results in terms of the physical parameters , 
mj 2 , Qi has been performed inE3. It yields qualitatively the same results. An 
additional source of deviation between the Feynman-diagrammatic and the RG 
result is the dependence of the diagrammatic result on the mass of the gluino, 
which does not appear in the RG result. Its variation gives rise to a shift of 
the diagrammatic result relative to the RG result of up to ±2 GeV. 
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